function parameter=parameter()
parameter.Rp=2*pi*10^3;  %pumping rate
parameter.s=[0;0;1];  %circular polarized pumping along z axis
parameter.T1=10^(-3)/(2*pi); %T1=1 ms
parameter.T2=2*10^(-5)/(2*pi); %T2=20 us
end
